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The reduced tension acd of the interface between the confined and the deconfined phase of SU (3) pure gauge 
theory is determined from numerical simulations of the first transfer matrix eigenvalues. At Tc = 1/Lt we find 
fed ~ 0.139(4)Te for Lt = 2. The interfaces show universal behavior because the deconfined-deconfined interfaces 
are completely wet by the confined phase. The critical exponents of complete wetting follow from the analytic 
interface solutions of a Z(3)-symmetric $ model in three dimensions. We find numerical evidence that the 
confined-deconfined interface is rough. 



The tension of the interface between the low 
temperature hadron phase and the high temper- 
ature quark-gluon plasma phase is an important 
parameter of the QCD phase transition. In the 
early universe the interface tension determines 
the nucleation rate of hadronic bubbles from the 
high temperature plasma. Thus, it sets the scale 
for spatial inhomogeneities in the baryon density, 
which may influence the primordial nucleosynthe- 
sis of light elements. To obtain an estimate for the 
interface tension in QCD we neglect the quarks 
and we concentrate on an SU (3) pure glue theory. 
Then the phase transition is first order, i.e. the 
confined and the deconfined phase coexist at the 
critical temperature Tc, and the interface tension 
is nonzero. It is convenient to define the reduced 
interface tension 



are three distinct high temperature phases, 
which coexist at temperatures T > Tc- Two 
different deconfined phases are separated by 
deconfined-deconfined interfaces with a reduced 
interface tension addiT). At Tc all four phases 
(three deconfined and one confined) coexist with 
each other. Thermodynamical stability requires 
<^dd{Tc) < 2crcrf- Frei and Patkos have sug- 
gested that 



(^dd{Tc) — 2a. 
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where F is the free energy of the interface and A 
is its area. 

In the pure glue theory the deconfinement 
phase transition is associated with the spon- 
taneous breakdown of the Z(3) center symme- 
try. Hence, as opposed to full QCD, there 



Then a deconfined-deconfined interface consists of 
two confined-deconfined interfaces with a macro- 
scopic layer of confined phase in between. One 
says that the deconfined-deconfined interfaces are 
completely wet by the confined phase. Com- 
plete wetting is a critical phenomenon of inter- 
faces (hence Tc is really a critical temperature 
although the bulk transition is first order). As 
one approaches the phase transition from above 
the thickness of the confined wetting layer 



zo oc (T - Tc) 
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diverges with a critical exponent tp, while the or- 
der parameter at the interface 

$i(0)oc (T-T,)^ (4) 

vanishes with another critical exponent (3. 
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In our case the order parameter is the Polyakov 
loop 

= Tr[-Pexp / dtAi{x,t)], (5) 
Jo 

where A^{x, t) is the anti-hermitean S'C/(3) gauge 
potential. The Polyakov loop ^{x) = + 
i^2{x) is a complex scalar field in three dimen- 
sions. To compute the values of the critical expo- 
nents one may construct an effective Z(3) sym- 
metric $^ theory with an action |^ 

= a|$|^ + - 3$2) + c|$|^ (6) 

The bare interfaces are solutions of the corre- 
sponding classical equations of motion. Fig.l 
shows the shapes $i(z) and $2(2^) of a planar 
deconfined-deconfined interface perpendicular to 
the z-direction deep in the deconfined phase (a), 
in a region where wetting sets in (b), and very 
close to the phase transition (c) where interface 
splitting and hence complete wetting is clearly 
visible. Close to the phase transition an ana- 
lytic interface solution yields zo « log(r— Tc) and 
$1(0) cx \/T — Tc. In particular, the thickness zq 
of the confined complete wetting layer diverges 
only logarithmically. The corresponding critical 
exponents are 

V'-O, /3 = i (7) 

Complete wetting is a peculiarity of the pure 
glue theory. In full QCD it does not occur be- 
cause the quarks break the Z(3) symmetry ex- 
plicitly, thus eliminating two of the three decon- 
fined phases. In numerical simulations of the pure 
glue theory wetting complicates the situation and 
makes the determination of the reduced interface 
tension acd more difficult. To avoid the compli- 
cations it is advantageous to impose C-periodic 
boundary conditions j^. Just like quarks they 
break the center symmetry and they reduce the 
number of deconfined phases from three to one 
0. With C-periodic boundary conditions the 
Polyakov loop is replaced by its charge conjugate 
when it is shifted a distance Li in the direction of 



Figure 1. The deconfined-deconfined interface. 
The solid line is <i>i and the dashed line is $2. 



the spatial unit vector 

$(f -Ki.e,) = = $*(f). (8) 

Hence, the two deconfined phases with non-real 
expectation values of the Polyakov loop are not 
consistent with C-periodic boundary conditions. 
Of course, one must use periodic boundary con- 
ditions in the euclidean time direction. 

To extract the value of acd we investigate the 
system in a cylindrical spatial volume with two 
short X- and y-directions of lengths Lx and Ly 
and a much longer z-direction. Close to Tc typi- 
cal configurations then consist of several blocks of 
confined and deconfined phases, aligned along the 
z-direction, and separated by confined-deconfined 
interfaces of area A = L^Ly spanned in the short 
X- and y-directions. The interface shape $i(z) 
(averaged over the short x- and y-directions) is 
shown in fig. 2 for a typical configuration on an 
8 X 8 X 128 X 2 lattice with C-periodic boundary 
conditions. The interfaces infiuence the spectrum 
of the transfer matrix in the z-direction. In par- 
ticular, the first two transfer matrix eigenvalues 
<o,i(z) = cxp{—Eq^iz) can be obtained from a 
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relation function of Polyakov loops separated in 
the z-direction. We have used lattices 4 x 4 x 64, 
4 X 6 X 64, 6 X 6 X 64, 6 X 8 X 96 and 8 X 8 X 128 
always with Lt = 2 lattice points in the eu- 
clidean time direction, and with several values of 
the Wilson couphng very close to Pc = 5.0933(7). 
Fig. 3a shows the energy Ei for simulations on a 
4 X 6 X 64 X 2 lattice with C-periodic boundary con- 
ditions. As expected the minimum of Ei occurs 
close to /3c, where fc — fd and hence a; = 0. Fig. 4 

Figure 2. A typical interface shape. 



2x2 matrix of transition amplitudes 

To lowest order of the dilute interface approxi- 
mation the amplitude for transitions between the 
confined and the dcconfincd phase is given by the 
one-interface contribution 

tcd{z) = / dzo 
Jo 

exp(— xzq)^ exp(— CTcrf^) exp(a;(z — zq)). (10) 

One integrates over the position zq of the in- 
terface. The Boltzmann factors exp(— a;2;o) and 
exp(a:(z — zq)) represent a block of confined phase 
of thickness zq and a block of deconfined phase 
of thickness z — zq, where x = \{fc ~ fd)A/T 
depends on the free energy densities fc and 
fd of the confined and the deconfined phase. 
The interface contribution 5 exp(— dcd^) is sup- 
pressed by the interface free energy F/T = cTcdA, 
and it also contains a factor S which describes 
capillary wave fluctuations of the interface. A 
rigid interface without capillary wave fluctua- 
tions has S — 1. Using similar expressions for 
the amplitudes tdd{z) and tcdz) and working 
to all orders of the dilute interface approxima- 
tion one finds an energy difference Ei — Eq — 
2^yx^ + 6'^ exp{—2acdA), which is minimal for 
a: = Putting £'o = one obtains at the 
minimum 

Ei^26exp{~acdA). (11) 

In the numerical simulations the energy difference 
is determined from the exponential decay of a cor- 



Figure 3. The spectra of the lowest energies. 



shows a fit of log Ei as a function of the interface 
area A — L^Ly. According to eq.(pl]) the slope 
of the curve determines the value of the reduced 
interface tension cr^d = 0.035(1) = 0.140(4)r2. 
The fluctuation factor S — 0.196(5) 1 indicates 
that the confined-deconfined interface is rough as 
it should be in the continuum limit. Note that 
rigid interfaces are lattice artifacts. 

When periodic boundary conditions are used 
the situation becomes more complicated, because 
then three deconfined phases are present and wet- 
ting may arise. In particular, there are four 
energies Eq, Ei, E2 and E3 which become de- 
generate at Tc in the infinite volume limit. In 
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Figure 4. The dependence of the energies on the 
interface area A. 



case of complete wetting (putting again Eq = 0) 
we find Ei^2 = \/ x'^ + 3(5^ cxp{—2acdA) — x and 
E:i = 2^a;2 +3(52exp(-2CTed^) (s). The energies 
El and E2 are degenerate because of Z(3) sym- 
metry and charge conjugation invariance. The 
energy E^ has its minimum at a; = where 



E3, = 2\/3(5exp(-(Tcd^) = 2Ei. 



(12) 



On the other hand, if wetting is incomplete (i.e. 
if addiTc) < 2(Tcd) the minimum of E3 occurs at 
X = — 7exp(— (Tdrf^) < 0, while the point where 
E3 = 2Ei^2 is shifted to x = 2^exp{—addA) > 
Q. The transfer matrix spectrum is then quali- 
tatively different. Fig. 3b shows the spectrum on 
a6x6x64x2 lattice with periodic boundary 
conditions. The minimum of E^ and the point 
where E3 — 2Ei,2 are only slightly shifted away 
from each other. In case of incomplete wetting 
one would expect a shift, but it would go exactly 
in the opposite direction. Therefore, our data are 
consistent with complete and not with incomplete 
wetting. In fact, using eq.(12) for an indepen- 
dent fit which is also shown in fig. 4, we obtain 
acd = 0.034(2) = 0.136(8)r2 and S = 0.202(6) in 
agreement with the C-periodic results. 

To summarize we have found evidence that the 
deconfined-deconfined interfaces are completely 
wet by the confined phase. Combining the C- 



periodic with the periodic results we obtain 

acd = 0.139(4)r2 (13) 

for the reduced confined-deconfined interface ten- 
sion. Our result is consistent with the one of the 
Boston group acd = 0.12(2)T^^ and larger than 
the one of the Helsinki group acd = 0.08(2)r2 0. 
As opposed to these methods we use the necessar- 
ily present finite size effects to extract the value 
of (Ted, instead of trying to avoid finite size ef- 
fects by going to larger lattices. In particular, a 
reliable extrapolation to the infinite volume limit 
is then possible, because we find agreement with 
our analytic finite size formulae. The fluctuation 
factor 5 — 0.199(4) indicates that the confined- 
deconfined interface is rough, as it should be in 
the continuum limit. Of course, our results are 
obtained with Lt — 2 and it would be interest- 
ing to get closer to the continuum limit, e.g. to 
Lt = 4. This is clearly feasible using our finite 
size technology, but it requires larger lattices like 
e.g. 16 X 16 X 256 x 4 because the phase tran- 
sition is then more weakly first order. Also high 
statistics is needed to be able to determine reli- 
ably the spectrum of the transfer matrix. The 
existing data already indicate that the interface 
tension is small at least in the pure glue theory. If 
this remains true in full QCD, only small spatial 
inhomogcncities are produced at the QCD phase 
transition and primordial nucleosynthesis is likely 
to take place under spatially homogeneous condi- 
tions. 
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